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ABSTRACT 

Spatially extended emission regions of Active Galactic Nuclei (AGN) respond to continuum variations, if 
such emission regions are powered by energy reprocessing of the continuum. The response from different parts 
of the reverberating region arrives at different times lagging behind the continuum variation. The lags can 
be used to map the geometry and kinematics of the emission region (i.e., reverberation mapping, RM). If the 
extended emission region is not spherically symmetric in configuration and velocity space, reverberation may 
produce astrometric offsets in the emission region photocenter as a function of time delay and velocity, de- 
tectable with future /ias to tens of /ias astrometry. Such astrometric responses provide independent constraints 
on the geometric and kinematic structure of the extended emission region, complementary to traditional rever- 
beration mapping. In addition, astrometric RM is more sensitive to infer the inclination of a flattened geometry 
and the rotation angle of the extended emission region. 

Subject headings: black hole physics — galaxies: active — quasars: general 



1. INTRODUCTION 

Reverberation mapping (RM) is a powerful tool to probe 
the structure of th e broad emission line region (BLR) in AGNs 
and quasars (e.g.. lBlandford & McKeelll 9821: (plitersonll 19931) 
without the need to resolve the BLR. The basic ideas are 
that the BLR is powered by photoionization by the (ioniz- 
ing) continuum from the accreting black hole and that it re- 
sponds to the variations of the continuum in a light crossing 
time. The response from different parts of the BLR will arrive 
at different time delays with respect to the continuum varia- 
tion. Thus by mapping a two-dimensional broad line response 
function in the time delay and velocity plane (velocity-d elay 
maps, e.g.. lBlandford & McKee1ll982l:lHorne et alj|2004 one 
can in principle recover the geometry and kinematics of the 
BLR. Over the past several decades, RM has proven to be 
a practical technique in studying the structure of BLRs. Al- 
though accurate velocity-time delay mapping is still lacking, 
RM studies have successfully measured average BLR sizes 
for several dozens of AGNs and quasars (e.g.. iKaspi et alj 
l2000l 120071: iPeterson et alHIOOl iBentz et alJl2009allbh. and 
in a few cases, crude ve locity-delay maps (e.g.. lDenney et alJ 
120091 IBentz eTani2010h . 

RM studies have shown that the typical BLR size R 
scales approximately with the AGN lu minosity L 05 . In th e 
latest version of the R — L relation (e.g.. IBentz et al.l l2009a). 
log(tf /light days) = -21.3 + 0.5191og(AL A (5100A)/erg S - 1 ). 
Thus for typical quasar luminosities (ALa(5100A) = 



10 45 ergs _1 , or bolometric luminosity Lboi ~ 10 46 ergs _1 ) at 
Z ~ 0.5, the BLR size is ~ 0.1 pc (~ 15 /ias). This scale is 3 
orders of magnitude smaller than the diffraction limit of 10m 
class optical telescopes (9 ~ tens of mas). Optical/near-IR 
interferometry with /ias resolution has yet to come. Thus 
reverberation mapping will remain one of the few practical 
methods to measure the BLR size, alon g with microlensing 
in gravitationall y lensed quasars (e.g., Morgan et al.l 120101: 
ISluse et alJfeOllI) . 

However, measuring the source photocenter positions can 
achieve a factor of ~ 1 / ^N p ^ oton enhance ment in precision 
compared wit h the image resolution (e.g.. lLindegren|[l978l: 
Bailevi n~998al) . where A'photon is the number of photons re- 



ceived in a bandpass. This means for 10 6 photons, the achiev- 
able astrometry precision is 3 orders of magnitude smaller 

than the image resolution; O ne a pplication of this idea is 

spectroastrometry (e.g., [Bailey 1998a), where differential as- 
trometric positions as a function of wavelength (veloci ty) can 
be used to probe otherwise unresolved sources (e.g., Bailey 
[19985 iGnerucci et al.ll201 ll) . 

When the astrometric precision approaches the angular 
BLR size, it becomes possible to study the BLR structure with 
astrometric signatures. A simple application would be to use 
spectroastrometry to place constraints on the BLR structure 1 . 
A more radical possibility, however, is to detect and model the 
wobble in the broad line emission photocenter due to reverber- 
ation to the continuum variations, since the different arrival 
times of the BLR response will cause shifts in the observed 
BLR photocenter. 

In this work we investigate the feasibility of combing the 
traditional intensity RM with astrometric information. Un- 
like spectroastrometry, the reverberation of the BLR will di- 
rectly induce shifts in the photocenter position of broad line 
emission as a function of time delay and velocity. As we will 
show below, the pattern of photocenter shifts is determined by 
the geometry and kinematics of the BLR, thus providing in- 
dependent constraints on the BLR structure, complementary 
to traditional intensity RM. Although our main focus is on 
the BLR, this method can be readily applied to the reverber- 
ation of the dust torus of AGNs (e.g., ISuganuma et alj|2006l) 
or other reverberation systems, where the required astromet- 
ric precision might be less stringent. We describe the basic 
formalism in $2] and demonstrate this method in §[3] and §|4] 
with simple models for the BLR. We briefly discuss the prac- 
tical issues with this method in $5] with an overview on the 
perspectives of achieving the required astrometric precision in 
the near future in 35.21 and conclude in §[6] 

2. BASIC FORMALISM 

Let us set up a 3-dimensional cartesian coordinate system 
(in the observer frame) centered on the continuum source 

1 For instance, in certain geometries of the BLR, the blue and red parts of 
the broad line emission may have offset photocenters. The detection of such 
photocenter offsets can put constraints on the BLR size and kinematics. 
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(e.g., the inner edge of the accretion disk, which we treat as 
a point source), and let the observer be at z = +00. The time 
delay of response from each point r = (x,y,z) in the BLR is 
then simply: 



r(r) = r-z = \/x 2 +y 2 +z 2 -z = r(l-cos0) , 



(1) 



where we use units that normalize the speed of light to c = 1 , r 
is the distance to the central source, and 9 is the angle between 
vector r and +z axis. Obviously the iso-delay surface (i.e., 
the surface with constant time delay) is a paraboloid for an 
observer located at infinity. Points on the +z axis, i.e., on the 
near side along the line-of-sight (hereafter los), have zero time 
delay with respect to the continuum variations. 

Assuming that the observed continuum increases at time t = 
by a constant amount f c for a certain amount of time Af , and 
goes back to normal at t = At, i.e., the change in continuum 
luminosity is 



Sc(t)=- 



t <0 

< t < At 
t > At . 



(2) 



The response in the broad line intensity is then a function of 
time t and los velocity v 



<W*,v)= / j(r)g(v,r,\r)6 c (t-T(r))drd\r , 



(3) 



where j(r) is the (assumed isotropic) responding volume 
emissivity of the emission region as a function of position, 
and 

g(v,r,w) = /(r,w)(5(v-w-n) (4) 

where (r,w) are the 3D coordinate and velocity vectors, 
/(r,w) is the normalized (e.g., J fdw = 1) 3D velocity dis- 
tribution at point r, determined by the kinematic structure of 
the BLR, and n = +z denotes the unit vector of the los. Note 
that Eqn. (01 is equivalent to the traditional RM equation in- 
voking the transfer function (i.e., the line intensity response 
for a ^-functio n continuum flare; Blandford & McKeel ll982l; 
lPetersonll993l) . However, for simplicity we shall avoid the us- 
age of transfer function and work directly in the time domain 
in terms of the observed line intensity and ph otocenter posi- 
ton fo r different velocity channels (see, e.g., iPancoast et al.l 
l20ll . 

If we further assume that the BLR is transparent to its own 
emission, i.e., neglecting absorption and scattering (as gener- 
ally assumed in RM), then the projected photocenter position 
of the responding emission with los velocity v as observed by 
the distance observer is 



xo(/,v) = 
Jo(?,v) = 



/ xj(r)g(v,r,w)6 c (t -r(r))drdw 
J j(r)g(v,r,w)S c (t-T(r))drdw 
J yj(r)g(v,r,w)6 c (t-T(r))drdw 
J j(r)g(v,r,w)S c (t-T(r))drdw 



(5) 



The projected photocenter position integrated over all los 
velocities is: 

Jxj(r)S c (t-T(r))dr J yj(r)S c (t - r(r)) dr 



J j(r)S c (t-r(r))dr 



Jj(r)S c (t-T(r))dr 



(6) 

In Eqn. (f2]), if cAt is much smaller compared with the typi- 
cal BLR size R, then the continuum light curve can be treated 
as a (square wave) 5-function. In this case the near side of the 



BLR completes reverberation (i.e., back to the steady level) 
before the far side reverberation reaches the observer. On 
the other hand, if cAt ^> R, then the continuum light curve 
is a step function, and the near side of the BLR will still be 
in the high state when the far side reverberation reaches the 
observer. In the following discussion, we will consider the 
two limiting cases where cAt <C R (a continuum flare) and 
cAt ^> R (a continuum jump). 

3. SIMPLE GEOMETRIC AND KINEMATIC MODELS FOR THE BLR 

To demonstrate the astrometric signals from RM, we use 
simple geometric and kinematic models for the BLR, and 
study the reverberation process. But our approach below is 
applicable to general cases. 

There are few observational constraints on the structure 
of the BLR. Spherical BLRs will not produce any photo- 
center offset. However, there is evidence that the BLR 
motion may be dominated by rotation, and the BLR may 
have a flattened geometry. Such evidence comes from the 
dependence of broad line width o n orientation angle in- 
ferred from radio pro perties (e.g., IWills & Brownd fl986t 
iJarvis & McLure 2006), or from mode ling the broad line pro- 
files (e.g.. iKollatschny & Zetzll 1201 lb . In the extreme case, 
the broad lines in < 10% AGNs show a characteristic double- 
peaked profile, commonly interpreted as arising from a disk 
geom etry with Keplerian rotation (e.g., lEracleous & Halpernl 
Il994h . On the other hand, radial motion of the BLR has 
been inferred from vel ocity-resolved rever beration mapping 
in several AGNs (e.g. IDennev et al.ll2.Q09b . and it has been 
suggested that the high-ionization broad line (such as Civ) 
may have a component originating from a disk wind (e.g., 
iMurrav et al.ll995t|Proga et al.l2000T) . It is possible that virial 
motion and radial motion coexist in the BLR. 

Motivated by these observations, we consider the follow- 
ing simple classes of geometry and kinematics for the BLR 2 . 
In all cases we assume constant density of the BLR, and we 
assume uniform reprocessing coefficient e across the entire 
BLR. Therefore the volume emissivity j(r) = eL c /(4irr 2 ) is 
proportional to 1/r 2 , i.e., spherical shells with the same thick- 
ness but different radii will have the same total reverberation 
intensity. Again, we emphasize that our approach is not lim- 
ited to these assumptions. We establish an intrinsic coordinate 
system (x'y'z') anchored to the BLR, which is rotated from the 
observer coordinate system (x,y,z) by Euler angles (a,/3,7). 

1. A triaxial ellipsoid configuration distributed within a 

, x' 2 y' 2 z' 2 , 
shell: r( < —r- + — 5- + — 5- < rr ut . We further consider 

a 1 b c 

a sub-class of this configuration, an oblate spheroid 
where a = b> c, which probably more represents a flat- 
tened BLR geometry. The velocity field is assumed to 
be Maxwellian, with ID dispersion equal to the virial 
velocity at each radius. Spherical shell geometries are 
also simplified versions of this model. 

2. A truncated circular razor- thin disk (R m < R < R out ) in 
the x'y' plane, with an inclination angle i from the los 
(z-axis in the observer frame). We assume a Keple- 
rian velocity field for the disk and the disk rotation is 
counter-clockwise viewed from z' = +00. 

2 Note that in this paper we do not intend to construct fully self-consistent 
models for the BLR, nor do we intend to reproduce the observed emission 
line profile with our BLR models. Rather, we use these idealized models to 
demonstrate the expected signals from reverberation mapping. More realistic 
BLR models can be easily implemented in our framework. 
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FIG. 1 . — Illustrative examples of BLR geometry and kinematics for the three non-spherical models described in ij3] and for the velocity-integrated reverberation 
amplitude 5 cra (0 and photocenter position Xo(t), The top panels show the cross section (regions filled with different patterns) in the xz. plane for three BLR 
geometries when the BLR intrinsic coordinate system x'y'z 1 is aligned with the observer coordinate system xyz. The red curves are typical parabolic iso-delay 
curves with constant r (in the same units as in Eqn.[T}. The bottom panels show the corresponding reverberation mapping signals. Left: an oblate spheroid model 
with a = b = 2c = 1 , r m = 0.5 and r out = 1 ■ Middle: a truncated disk model with R m = 0.5 and R OM = 1 . Right: a biconical outflow with r m = 0. 1 , r 011 t = 1 .0, and 
an opening solid angle Q = tt/2. In all examples, the intrinsic coordinates of the BLR (x'y'z 1 ) are rotated about the y-axis in the observer frame by an inclination 
angle i, such that the projected BLR in the xy plane has planar symmetry about the .v-axis (hence yo(t) = 0). For each model we compute for three inclination 
angles i = 0° .45° ,60°, and two continuum variation timescales At — 0.2 and 0.3. In the (' = 0° cases, there is no photocenter offset as a function of time due to 
symmetry. We define the maximum-possible reverberation intensity as the one where the entire BLR is reverberating. 

kinematics to sufficient numeric accuracy 3 compared with in- 
tegrating Eqns. (O and (0. For simplicity we adopt the units 
c=G = M = 1, where G is the gravitational constant and M is 
the central black hole mass. 

4. RESULTS 

In the simplest case of a thin spherical shell (with radius 
R) model for the emission line region, the time evolution of 
the reverberation in tensity has simple analytical forms (e.g., 
iBahcall et al] 1 19721) . which we verified with our numerical 
approach. However, due to the spherical geometry, the photo- 
center does not show any offset as a function of time. 

Fig. [T] shows several examples for the three non-spherical 
models described in $3] f° r the velocity-integrated case. Since 
all three models have rotational symmetry, we need only focus 

3 In all the calculations shown below we have performed convergence tests 
to make sure that the grid is fine enough to capture the results to < 1 % relative 
accuracy. 



A biconical configuration, which is rotationally sym- 
metric about the z'-axis with an opening solid angle f2, 
and a radial extent from r m to r out . The bicone has an 
inclination angle i from the los (z-axis in the observer 
frame). We assume a biconical radial outflowing ve- 
locity structure. The detailed outflow velocity structure 
is not important in this work, as we will only consider 
the blueshifted and redshifted halves of the broad line 
profile. 



To study reverberation mapping in emission line intensity 
and photocenter offsets for arbitrary geometry and kinemat- 
ics, we discretize the emission region with a fine cartesian 
grid of cells (x,-,y,,z,) and assign velocities to each cell. At 
each time t , we determine the cells that are reverberating, i.e., 
$c(t — T~(xi,yi,Zi)) 7^0, and compute the reverberation quanti- 
ties 6 em , xo and yo in Eqns. (f3]l and <(5j with direct summation. 
This discrete method can efficiently handle any geometry and 
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FIG. 2. — Velocity-resolved evolution of reverberation intensity and photo- 
center positions, for the same disk model as in Fig.[T]with a los angle (' = 45° 
and two continuum variation timescales At = 0.2, 3. Due to symmetry, the 
intensity <5 cm and .v-positions of the photocenter have identical time evolu- 
tion for the blueshifted and redshifted halves of the disk, which are denoted 
by the black lines. On the other hand, the y photocenter positions have time 
evolutions denoted by the blue and (thicker) red lines for the blueshifted and 
redshifted halves, which are symmetric about the ;t-axis. 



FIG. 3. — Velocity-resolved evolution of reverberation intensity and photo- 
center position xo(t), for the same bicone model as in Fig.[T]with a los angle 
(' = 45° and two continuum variation timescales Af = 0.2, 3. The blueshifted 
cone (denoted by blue lines) and the redshifted cone (denoted by thicker red 
lines) have quite different time evolution in reverberation intensity <5 cm and 
photocenter position xo(t), reflecting the geometry and kinematics of the bi- 
cone BLR structure. The photocenter position yo(t) = due to symmetry in 
the projected xy plane. 



on the photocenter offset in the x-direction, i.e., the rotation 
of our model BLR preserves symmetry about the x-axis in 
the projected xy plane. Specifically, the intrinsic coordinate 
system of the BLR (x'y'z 1 ) is rotated about the y-axis in the 
observer frame by inclination angle i. 

For each model we compute for three orientations i = 
0,45°, 60° and two continuum variability durations Af = 0.2 
and 3. In all cases the maximum BLR radius is set to be 1. 
It is clear from Fig. Q]that similar information is encoded in 
the photocenter reverberation maps as in the traditional inten- 
sity reverberation map. However, the photocenter reverber- 
ation maps are more sensitive to the inclination of the flat- 
tened BLR geometry than the intensity reverberation maps. 
The maximum photocenter offset is on the same order as the 
flux-weighted average size of the BLR. We will return to these 
points in $5] 

In the examples shown in Fig. [T] a generic feature of the 
photocenter reverberation is the transition from the near-side 
photocenter to the far-side photocenter. If the duration of the 
continuum variation Af is short compared with the BLR size 
(e.g., a flash), this transition is continuous. On the other hand, 
if the duration Af is long, then there will be a period with no 
photocenter offset (i.e., the entire BLR is lit up). The spheroid 
geometry and the disk geometry result in similar astrometric 
RM signatures. This is expected because the spheroid model 
used is close to a flattened disk geometry. In both cases we 



see a gradual rise followed by a sharp cutoff around the max- 
imum photocenter offset after the continuum returns to the 
base value. However, in the bicone geometry, the sharp cutoff 
occurs later past the peak photocenter offset. This difference 
is caused by the different edge geometry in these models - in 
the disk and spheroid cases, fewer and fewer materials are lo- 
cated at the far side of the emission region, which is opposite 
in the bicone case. 

Another point worth noting is that for inclined BLRs, the 
total line intensity reverberation 8 em (t) looks similar in dif- 
feren t models, thus doe s not offer much differentiating power 
(e.g.. lHorne et al.|[2004h . However, the different behaviors in 
photocenter reverberation xo(f) for the flattened geometry and 
for the bicone geometry might be useful in distinguishing the 
two geometries. 

The astrometric reverberation mapping signals vanish when 
the iso-delay surface is axisymmetric about the z-axis. For 
instance, an edge-on disk will not produce photocenter offset 
as a function of time delay for the velocity-integrated broad 
line flux. Fortunately, an astrometric signal will still arise if 
we have velocity-resolved reverberation mapping: the blue 
side of the reverberating broad line has a photocenter shifted 
from that of the red side. Thus the astrometric feature can 
determine the angular momentum direction of the disk-like 
BLR, while the traditional intensity reverberation mapping is 
unable to derive such information. 



Astrometric RM 



5 



Similarly, velocity-resolved reverberation mapping in flux 
amplitude and in photocenter offset is a powerful tool to fur- 
ther constrain the velocity field of the BLR. Below we use the 
three models to demonstrate the different reverberation behav- 
iors in the blueshifted and redshifted halves of the line profile. 

1. Spheroid models. Since we assumed a random 
Maxwellian velocity distribution, the red part and blue 
part of the emission line profile have identical behavior 
in reverberation signals. 

2. Disk models. We consider the disk example of i = 45° 
shown in Fig.Q] In this case the redshifted part (vi os < 
0) and blueshifted part (vi os > 0) of the disk will be- 
have exactly the same in terms of S em (t) and xo(t), ex- 
cept that the amplitude of S em (t) is reduced by half. 
However, since we are now looking at half of the disk, 
the y photocenter position of the redshifted/blueshifted 
disk changes as a function of time, as indicated by the 
red/blue lines in Fig. [2] 

3. Bicone outflow models. We consider the bicone exam- 
ple of i = 45° and ft = tt/2 shown in Fig. [T] In this ex- 
ample the near-side cone has blueshifted velocity while 
the far-side cone has redshifted velocity 4 . We plot the 
results for both cones in Fig. [3] The line intensity S em is 
already different for the blueshifted half and redshifted 
half, since the two cones have distinct los velocities: we 
see intensity reverberation from the blueshifted cone 
(near-side) earlier than the redshifted cone (far-side). 
The time evolution of the x photocenter position is also 
different for the blueshifted and redshifted cones. The y 
photocenter positions for the two cones are always zero 
due to the symmetry of this example. 

In the traditional case, velocity-resolved intensity RM can 
effectively distinguish outflow/inflow kinematics from virial 
motion (either in a Keplerian disk or in randomly orientated 
orbits) of the emission line region. The distinction between 
the two kinds of kinematics is also clearly reflected in the 
velocity-resolved astrometric RM, i.e., the time evolution of 
the blueshifted and redshifted photocenters is asymmetric in 
the outflow/inflow case, and symmetric in the case of virial 
motion (cf. Figs.[2]and[3]l 

5. DISCUSSION 

5.1. Practical Concerns with this Method 

The representative examples shown in §|4] indicate that the 
expected astrometric signals are on the order of the BLR size. 
Then for typical quasar luminosities of L\, \ ^ 10 46 ergs _I , 
we expect such astrometric signals to be on the order of 
tens of [ias at z ~ 0.5. For more luminous objects and/or at 
lower redshifts, the expected astrometric signals will be larger. 
For example, the quasar 3C 273 (z = 0.158, AL A (5100A) = 
8.6 x 10 45 ergs _1 ,V ~ 12.8) has a meas ured H/3 BLR size 
of ~ 300 light days (Be ntz et"aT1 l2009ab . corresponding to 
~ 90 /Lias at its redshift. If we further consider reverberation of 
the dust torus, then the expected astrometric signals would be 
even larger, e.g., > 100 /ias for quasar luminosities. The sizes 
of the BLR and torus approximately scale with L° 5 , thus for 

4 If the half-opening angle of the bicone is larger than 90° — i then the 
near/far side of the bicone will also contribute to redshifted/blueshifted ve- 
locities. 



Seyfert luminosities, /ias astrometric precision might be nec- 
essary for astrometric RM. 

However, even if the astrometry requirement is fulfilled (see 
{35.21 i. there are still several practical issues to consider. These 
will affect the target selection and observing strategy for a 
successful monitoring program. Detailed characterization of 
the adverse effects and methods to mitigate them are beyond 
the scope of this paper, and we only give a brief discussion 
here as a general guideline for future work. 

First, the monitoring should be carried out in optical/near- 
IR to observe the restframe UV to near-IR broad emission 
lines (e.g., Civ, Mg II, Balmer and Paschen lines) and adja- 
cent continua. Tunable narrow-band filters are desired for the 
broad line imaging and reliable continuum subtraction. High 
photometric accuracy is generally required, and to detect and 
model the astrometric signals we need good cadence of mon- 
itoring data. The monitoring time baseline could span weeks 
to years, depending on the target. Observations of repeated 
reverberation events can be combined to provide improved 
model constraints, as long as the monitoring time span is 
much shorter than the dynamical timescale (rd yn ~ R / Vfwhm) 
of the BLR. 

Measuring the photocenter of the reverberating part of the 
BLR to the nominal astrometric accuracy poses a more dif- 
ficult challenge. It requires the subtraction of the continuum 
and the steady-state BLR emission and subsequent measure- 
ment of the photocenter of the variable BLR emission without 
degrading the astrometric accuracy much. The use of narrow- 
band filters helps with this procedure, but also increases the 
integration time to achieve the required photometric precision. 
Take the broad H/3 line for example, the typical restframe 
equivalent width is ~ 70 A (e.g. JShen et al.ll2~01 ll) . Thus for a 
~ 100 A narrow -band filter centered on the redshifted H/3 line 
at z ~ 0.5, the enclosed H/3 line flux is ~ 50% of the total flux. 
The host galaxy contribution in this bandpass is negligible for 
quasar luminosities AL A (5100A) > 10 45 ergs _1 , and reaches 
~ 20 % of the total con tinuum at AL A (5100A) = lO^ergs" 1 
(e.g. JShen et"aT][20Tlh 5 . With SNR> 100 for the PSF core, 
> 10% fractional changes in the broad line flux can be de- 
tected. However, the centering accuracy of the variable BLR 
emission degrades as cr me as,var ~ cr me as//var, where CT me as is the 
nominal centering accuracy for the target and / var is the frac- 
tion of the variable flux (from the reverberating BLR) to the 
total flux 6 . Thus larger amplitudes in the continuum vari- 
ations (hence larger reverberation amplitude) are always fa- 
vored over smaller amplitudes, just as for traditional intensity 
RM. Optimal image subtraction methods may be required to 
isolate the v ariable BLR emission within the bandpass (e.g., 
lAlardll2000b . 

Finally, the real situation may deviate from the idealized 
case assumed in reverberation mapping, which will inevitably 
lead to complications. With better and better RM data sets, 
however, it is possible to incorporate additional physical in- 
gredients, such as photoionization processes, absorption and 
scattering, more complex geometry/kinematics, etc., in the 
framework of combined intensity and astrometric RM. 

5.2. Precision Astrometry in the Optical/Near-IR 

5 These host contamination estimates were based on SDSS spectra with 
3" diameter fibers. For our purposes, the extraction aperture would be much 
smaller, so host contamination is expected to be significantly less. 

6 For bright targets, however, the astrometric accurac y wi ll be limited by 
systematic errors rather than by centering accuracy (see j|5.2t . 
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We now discuss realistic timescales for achieving the as- 
trometric precision (/ias to tens of /ias) required to detect the 
astrometric reverberation mapping signals. Such single-epoch 
astrometric precision is already under serious considerations 
for ground-based and space-based facilities. 

5.2.1. ground-based facilities 

The recent development of AO-assisted, high-angular reso- 
lution, near-infrared observations with large, single-aperture 
telescopes has dramatically improved the astrometry preci- 
sion. Single-measurement precision of ~ 200 -300 /ias has 
been routinely achieved on 8- 10m cl ass telescopes (e.g., 
lLazorenko et~aT]|2009t IFritz et all 1201 Oh . < 100 /ias single- 
epoch precision (stable over 2 months) has been demonstrated 
for th e Palomar 200 inch H ale telescope in 2-minute expo- 
sures ( ICameron et al.ll2009h . for stars with K s < 13. At this 
level of precision, there are several important instrumental, 
atmospheric and astrophysical effects that need to be taken 
into account; and in most cases, optimized procedures or al- 
gorithms must be undertaken to minimize these effects ( e.g., 
ICameron et ai1l2009l: IFritz et alJl20p ITrippe et alJl20Toh . 

In the study by ICameron et all (2009)," it was shown that 
the dominant systematic error is due to atmospheric differen- 
tial tilt jitter. This is the stochastic and achromatic fluctuation 
in the relative displacement of the target and the reference 
star due to the different columns of atmosph eric turbulence 
travers ed by the light from the two objects. ICameron et alJ 
(2009) demonstrated that the tilt jitter effect can be largely 
mitigated by an optimal weighting scheme that uses a grid of 
reference stars to measure the pair-wise distances from the tar- 
get. This optimal estimation technique can essentially reduce 
the error due to tilt jitter, ajj, below the centering accuracy 
Cmeas, which is mostly determin ed by aperture size, photon 
noise, and properties of the PSF. ICameron et all (12009b also 
provided an empirical formula to predict the astrometric per- 
formance of a single-conjugate AO system: 



2 _ 2 
tot meas 



+ Otj = 



1 .4 sec 



- mM N) Id 



2 mas 



N 



5m 



(7) 



where t is the integration time, D is the telescope aper- 
ture diameter and N is the number of reference stars. This 
formula suggests that for 20-minute exposures, the single- 
measurement precision can achieve ~ 20 //as for lOm-class 
telescopes, and < 10 /ias for 30m-class telescopes, with only 
a few reference stars in a 25" x25" FOV. Note that with these 
large-aperture telescopes and a reasonable amount of expo- 
sure time (e.g., 30-60 min), the magnitude limit achievable 
with the above astrometri c precision can be as faint as K s ~ 
18 by extrapolation (e.g., ICameron et alj|2009l : iTrippe et ail 
2010), relevant for reverberation monitoring of AGN contin- 
uum and broad line fluxes. 

The above estimation might be optimistic. The actual astro- 
metric performance will inevitably depend on t he telescope 
and A O systems, as well as target properties. Trip pe etall 
(2010) performed a detailed study on the limiting factors on 
the astrometry precision for the Multi-adaptive optics Imag- 
ing CAmera for Deep Observations (MICADO) proposed for 
the 42-m European Extremely Large Telescope (E-ELT). The 
dominant error terms include: instrumental geometric dis- 



tortion (Tdist ~ 30 /ias, chromatic atmospheric differential re- 
fraction (CDR) (Tcdr ~ 20 /ias, atmospheric differential tilt 
jitter otj ~ 10 /ias (with few minutes of integration), aniso- 
planatism of the multi-conjugate AO system cr an i so ~ 8 /ias, 
using galaxies as astrometric references in high galactic- 
latitude fields (Tg a i ax i es ~ 20 /ias, and calibration of the pro- 
jected pixel scale g sca ie 5 , 10 /ias. Combining these error 
terms, ITrippe et all d2010h estimated a total error budget of 
c.sys ~ 40 /ias. In practice, some of these error terms can be 
reduced. For instance, ocdr can be effectively suppressed 
if narrow-band filters are used or special algorithms and ob- 
serving strategies are deployed to correct for the CDR (e.g., 
ICameron et alj|2"0"09h . Note that for BLR reverberation map- 
ping, it is preferable to use narrow-band filters to cover the 
emission lines and nearby continuum. o- ga i ax i es can also be re- 
duced for AGN/quasar fields where there are a suitable num- 
ber of reference stars. It is also possible to correct instru- 
mental distortion by dedicated calibration proced ures to con- 
trol d dist down to the level of ~ 10-30/ias (e.g., Trip pe et"aT] 
2010). Thus the range of single-epoch astrometry precision 
for this particular instrument on E-ELT is ~ 20 - 40 /ias, al- 
though the lower bound would be very challenging to achieve. 
We expect similar performance for other 30m-class tele- 
scopes. 

An alternative route to single-aperture precision astrometry 
is long-baseline astrometry that operates in the very-narrow- 
angle regime (angular separation < 30"), which uses interfer- 
ometry in optical/ near-IR to achieve hig h relative astromet- 
ric precision (e.g., Shao & Colavita1ll992h . For instance, the 
Very Large Telescope Interferometer ( VLTI) instrument is de- 
signe d to deliver an astrom etric precision at the ~ 30 -40 /ias 
level dvan Belle et al. Il2008l) . In the short term, however, the 
astrometric precision is likely limited, by the accuracy in de- 
termining the in terferometer baseline, to ~ 20-30/ias (e.g., 
Ivan Belle! l2009h . Moreover, current implementation of op- 
tical to near-IR interfero metry usually require s bright target 
magnitude K < 13 (e.g.. Ivan Belle et alll2008h . hence is not 
yet suitable for AGN monitoring except for the brightest ob- 
jects (such as 3C 273). 

To summarize, ground-based facilities are promising to de- 
liver < tens of /ias single-epoch astrometry precision within 
a decade or two, and in the most optimal cases, ~ 10/ias pre- 
cision. Single-aperture astrometry with large telescopes and 
AO (currently operating in the near-IR, with potential exten- 
sion to the optical) is probably better suited for AGN targets 
than long-baseline interferometry at this stage. This preci- 
sion can be achieved for relatively short integration (< 1 hr) 
and faint target limit K s < 1 8 relevant for AGN reverberation 
mapping. But to achieve this goal, dedicated calibration algo- 
rithms and observing strategies, along with systematic control 
of the telescope and instrument system, will all be necessary. 
This level of astrometry precision is suitable for BLR rever- 
beration mapping of nearby AGNs and bright quasars, and 
sufficient for AGN dust reverberation mapping. 

5.2.2. space-based facilities 

The Hubble Space Telescope (HST) is currently the only 
space-based facility that can perform astrometry at the < 
1 mas precision (e.g.. I Anderson & King 120001 iBenedict et alJ 
2003). The future James Webb Space Tele scope (JWST) can i n 
principle deliver a similar precision (e.g., Diaz-Miller 20071) . 
Thus both HST and JWST are less favorable for precision 
astrometry compared with ground-based facilities discussed 
above. 
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GAIA dPerrvmanet all 1200 ll) is a space-based astrometry 
mission scheduled to launch in 2013, which can deliver 150 
(400) /ias single-measurement astrometric precision for V < 
16 (18) targets. This precision is not adequate for astrometric 
reverberation mapping of the BLRs, and barely enough for 
AGN dust torus reverberation. 

Perhaps the best facility for astrometric AGN reverberation 
mapping is a mission similar to the SIM/SIM Lite mission 
dUnwin et alj|2008l) . which is capable of achieving ~ 10/ias 
single-measurement precision for V < 18 targets with a space- 
based 6-m baseline optical interferometer 7 . In fact probing 
the <pc spatial extent of AGN/quasar jets and BLRs is one 
of the science cases for this mission, using the same tech- 
nique as spectroastrometry. Astrometric reverberation map- 
ping would be an ideal application with similar facilities for 
pointed multi-visit observations, yet with much more con- 
straining power on the BLR (and torus) geometry and kine- 
matics than spectroastrometry. Although this mission was 
canceled for the current decade, t he technolo gy for such a 
mission is already in place (e.g.. IShaol l2010h . In light of 
the exciting science that will be enabled with /ias astrome- 
try (especially with the driving force for exoplanet science), it 
is reasonable to anticipate similar missions to be approved in 
the next decade. 

6. CONCLUSIONS 

In this paper we outlined a simple idea of astrometric rever- 
beration mapping: observing the astrometric "wobble" of the 
BLR photocenter in reverberation events. We demonstrated 
its potential in constraining the geometry and kinematic struc- 
ture of the BLR with simple BLR models. 

The virtue of astrometric RM is that it provides an indepen- 
dent set of time series that can be combined with the intensity 

7 http://sim.jpl.nasa.gov/index.cfm 

8 Other methods to in fer the inclina tion of a disk geometry of the BLR 
include polarimetry (e.g., Li et al. 2009) and double-peaked broad line profile 
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RM data to constrain the BLR model. In addition, it is more 
sensitive to the inclination of a disk-like geometry 8 and infer- 
ence of its angular momentum direction, than intensity RM 
alone. Although the BLR is still unresolved, the spatial infor- 
mation from precision astrometry for the reverberating parts 
of the BLR greatly facilitates using reverberation mapping to 
probe the BLR structure. Once we have a set of high qual- 
ity continuum and broad line flux monitoring time series with 
precision photometry and astrometry, we can fit the data in 
the time domain using model intensity light curves and pho- 
tocenter evolution. General Markov C hain Monte Carlo ap- 
proaches with Baye sian inference (e.g., iPancoast et alJl2l)llt : 
iBrewer et alJ 1201 lh can be applied to search the parameter 
space of different model families and to quantify model un- 
certainties. 

While conceptually simple, this method is not yet ready 
for immediate applications. We discussed realistic timescales 
for achieving the required astrometric precision (/ias to tens 
of /xas), and concluded that ground-based large-aperture tele- 
scopes with AO offer the best hope to apply this method to 
luminous quasars and dust torus RM within a decade or two, 
and a space-based interferometer (a SIM/SIM-Lite like mis- 
sion) with /ias astrometric precision in the next decade or two 
may extend this method to the Seyfert regime. 
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